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ABSTRACT 

We have explored the redshift evolution of the luminosity function of red and 
blue galaxies up to 2; = 3.5. This was possible joining a deep I band compos- 
ite galaxy sample, which includes the spectroscopic K20 sample and the HDFs 
samples, with the deep Hab = 26 and Kab = 25 samples derived from the deep 
NIR images of the Hubble Deep Fields North and South, respectively. About 
30% of the sample has spectroscopic redshifts and the remaining fraction well- 
calibrated photometric redshifts. This allowed to select and measure galaxies in 
the rest-frame blue magnitude up to z ~ 3 and to derive the redshift evolution of 
the B-band luminosity function of galaxies separated by their rest-frame U — V 
color or specific (i.e. per unit mass) star-formation rate. The class separation 
was derived from passive evolutionary tracks or from their observed bimodal dis- 
tributions. Both distributions appear bimodal at least up to 2; ~ 2 and the locus 
of red/early galaxies is clearly identified up to these high redshifts. Both lumi- 
nosity and density evolutions are needed to describe the cosmological behaviour 
of the red/early and blue/late populations. The density evolution is greater for 
the early population with a decrease by one order of magnitude at 2; ~ 2 — 3 with 
respect to the value at z ~ 0.4. The luminosity densities of the early and late 
type galaxies with Mb < —20.6 appear to have a bifurcation at z > 1. Indeed 
while star-forming galaxies slightly increase or keep constant their luminosity 
density, " early" galaxies decrease in their luminosity density by a factor ~ 5 — 6 
from z ~ 0.4 to z ~ 2.5 — 3. A comparison with one of the latest versions of the 
hierarchical CDM models shows a broad agreement with the observed number 
and luminosity density evolutions of both populations. 
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1. Introduction 

The redshift evolution of the galaxy luminosity function in different bands is a powerful 
probe of the main physical processes governing the assembly and subsequent activity of 
galaxies. In particular, the analysis of the luminosity function of galaxies having different 
spectral and/or morphological properties can help to disentangle different histories of star- 
formation. 

In the emerging picture, galaxies form hierarchically through mergers of pre-existing 
smaller galaxies in the context of Cold Dark Matter (CDM) scenario. In this picture early 
type galaxies, which are in general among the brightest and massive galaxies in the universe, 
tend to be the result of relatively late assembly. This appears in contrast with the old age 
inferred from their red colors and with recent observations by Cimatti et al. (2004) and 
Glazebrook et al. (2004) showing the presence of some massive galaxies at high redshift. 

For this reason it is very important to analyze the abundance as a function of redshift of 
early type galaxies selected on the basis of their morphological appearance or spectral energy 
distributions (SED). This is a difficult task both in the observational and theoretical context: 
on the one side we need large and deep surveys of NIR selected galaxies to sample the rest- 
frame visual band at least up to 2; ~ 3; on the other side we need sophisticated MonteCarlo 
or hydrodynamical codes for galaxy formation and evolution in the CDM scenario. 

Most of the information on the statistical analysis of the red/early or blue/late type 
galaxies comes from local samples (e.g. Blanton et al. 2003, Baldry et al. 2003, Brinchmann 
et al. 2004). In these studies a clear bimodality is appearing in the color distribution of 
galaxies or in their specific (i.e. per unit mass) star-formation rate. 

This bimodality separates the locus of the star-forming galaxies from the early popu- 
lation (early type spirals, SO, EO) and thus represents a useful tool to study the redshift 
evolution of galaxies having colors consistent with those of the early types independently of 
any morphological selection. In the local universe these different definitions appear to be 
consistent with each other (see e.g. Schweitzer and Seitzer 1992). 

In general it is difficult to extract information about the star formation histories of 
individual early type galaxies since the stellar population evolving passively can become red 
in a time much shorter than the age of the galaxy. A complementary approach is given by 
the analysis of the redshift evolution of the luminosity function of galaxies having different 
colors or specific star-formation rates. 

Recent attempts in this respect have been made on the local Sloan sample e.g. by Hogg 
et al. (2002), Blanton et al. (2003), Baldry et al. (2004) and at intermediate redshifts 
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by Bell et al. (2004). Preliminary results show that the rest-frame U — V distribution of 
galaxies is bimodal up to 2; ~ 1 even for the field population. This allows a separation of the 
galaxy population in red/early and blue/late galaxies. In particular, the rest-frame B band 
luminosity function of the red galaxies shows a rather flat slope and evolves only mildly in 
qualitatively good agreement with the ACDM scenario. 

An attempt to follow at higher redshifts the evolution of galaxies of different spectral 
types has shown in particular that the number density of early type galaxies decreases with 
increasing rcdshift (Wolf et al. 2003). However the authors used fixed spectral energy 
distributions to separate the galaxy spectral types neglecting evolutionary effects due to the 
passive stellar evolution with redshift. Moreover, the lack of near infrared imaging in their 
sample has limited the analysis to 2; < 1. 

The extension of this information at high redshifts requires indeed deep NIR surveys 
to get an estimate of the rest-frame U — V color. The J,H band images in the HDFN 
and the J,K band images in the HDFS are among the deepest NIR images present in fields 
where extensive multicolor optical imaging is already available and have been used to derive 
accurate photometric redshifts (Dickinson et al. 2003, Labbe et al. 2003, Fontana et al. 
2003, Poh et al. 2003). 

In this context, we have already used in Poll et al. (2003, paper I) a composite sample 
of NIR selected galaxies that, although relatively small, has unique characteristics in many 
respects. It is a complete sample selected in the rest-frame B band from low to high redshifts, 
just thanks to the deep NIR HDF south images. It includes the area of the K20 survey 
(Cimatti et al. 2002) with a large fraction of spectroscopic redshifts at the bright end and 
well tested photometric redshifts at the faint end. This composite sample has been used to 
derive the evolution of the global rest-frame B-band luminosity function up to 2; ~ 3. The 
main result was the presence of appreciable luminosity evolution of the bright end of the LF. 
Such a brightening could be produced by enhanced star-formation at high redshifts triggered 
by galaxy encounters in a CDM scenario (see Menci et al. 2004, Nagamine et al. 2004). 
This kind of evolution is dominated by the star-forming population, which is abundant at 
high redshift. 

On the other hand, we have also explored the evolution in mass of the sample, exploring 
in particular the evolution of the mass function up to intermediate redshifts and the mass 
density up to 2; ~ 3 (Fontana et al. 2003, 2004). This approach allowed us to follow the 
cosmological evolution of galaxies in a way less biased by episodes of starburst activity. At 
^ < 1 little evolution is observed in the galaxy mass function derived from the K20 sample, 
while at higher redshifts a more pronounced decrease in the stellar mass density was found 
in the HDF-S up to ^ ~ 3 in broad agreement with ACDM models (see also Rudnick et al. 
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2003). However, the derived fraction of mass density with respect to the local value at these 
high redshifts shows field-to-field variations by about a factor 2 if we compare the results 
obtained in the HDF-S field with those found in the HDF-North by Dickinson et al. (2003). 

To explore how different the star formation history can be in different galaxy popula- 
tions, we derive in this paper the number and luminosity density evolutions of the galaxies 
as a function of their color or specific (i.e. per unit stellar mass) SFR (SSFR) up to 2; ~ 3, 
and we compare the density evolution of the red/early or blue/late populations with that 
predicted by a MonteCarlo rendition of the Menci et al. (2004) ACDM model. In the paper 
all the magnitudes arc in the AB system and an Q\ = 0.7, Qm = 0.3, and Hq = 70 km s~^ 
Mpc~^ cosmology was adopted. 

2. The galctxy sample 

Our composite sample of galaxies covers, with similar depth, a range spanning from 
the UV down to the optical/NlR wavelengths and is described in detail in paper 1. Here 
wc summarize its general characteristics and main improvements. The galaxies arc selected 
in the / or NIR J, K bands depending on the redshift interval as explained below. A first 
photometric catalog has been derived combining the optical images in the HDF-S (Casertano 
et al. 2000) with the deep ESO/VLT NIR images obtained in the framework of the FIRES 
project (Labbe et al. 2003). Our infrared catalog reaches formal 5 sigma limits in an aperture 
of 1.2" at Jab = 25.3, Hab = 24, and Kab = 25. In the I band the HDF-S galaxies have 
been selected down to Iab = 26.2. In the HDF North we used the same Kab < 23.5 catalog 
as in paper I. In addition to the sample used in paper I, we adopted the HDF-N multicolor 
catalog of galaxies selected in the NICMOS H band by Dickinson et al. (2000) which is 
highly (> 95%) complete down to Hab < 26 or Iab < 25.8, where the relatively bright 
hmit in the I band was required to ensure detection in the NIR bands. We used also a 
brighter sample of galaxies selected in two fields centered around the QSO 0055-269 and in 
the Chandra Deep Field South (Giacconi et al. 2001) for a total of ~ 68 arcmiv?. This 
sample was used to select the targets for the so called K20 {Kab < 21.9) spectroscopic 
survey (Cimatti et al. 2002). The catalogs in the latter two fields have been used to derive 
complete (> 95%) samples down to Iab = 24.5, 25, respectively. Detection in the NIR bands 
is ensured for all the galaxies in the composite sample. 

To obtain the total fluxes of the objects in the chosen reference band, SExtractor Kron 
magnitudes (Bertin & Arnouts 1996) have been adopted for bright sources. For faint sources 
in ground-based images we estimated aperture magnitudes with aperture corrections com- 
puted estimating the flux losses outside the adopted aperture with respect to the Kron 
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magnitudes for relatively bright sources. The reliability of the correction adopted has been 
tested by means of simulations as performed in paper I and in Poli et al. (2001). No appre- 
ciable systematic losses of the total flux for the faintest sources were found in addition to 
the well known losses by 5-10% typical of the Kron magnitudes. Moreover, Bouwens et al. 
(2004) have recently compared images obtained with HST/ACS in the Ultra Deep Field with 
those obtained in the same area from the shallower ACS-GOODS survey. The main result is 
that there are no specific biases due to surface brightness dimming in the evaluation of the 
typical galaxy sizes and consequently on the estimate of their total magnitudes up to high 
redshifts. To cope with the inhomogeneous quality of the different images, we produced a 
set of seeing— corrected color magnitudes that have been then scaled to total magnitudes in 
all the bands (see e.g. Vanzella et al. 2001). This ensures that colors - used to estimate the 
redshifts - are measured in the same physical region at all wavelengths, independently of the 
seeing FWHM present in each image. A fraction of 29% of the sample has spectroscopic red- 
shifts, mainly from the K20 survey (Cimatti et al. 2002) and from the spectroscopic follow 
up in the GOODS/Chandra field (Vanzella et al. 2004), and in the HDF north (e.g. Cohen 
et al. 2000). For all the other galaxies, rehable photometric redshifts have been derived with 
a homogeneous technique, described in Fontana et al. (2000) and in Cimatti et al (2002). 
The model adopts several galaxy ages and exponentially declining star formation histories, 
and includes Lyman absorption from the intcrgalactic medium and diist absorption with 
different amounts up to E{B — V) — 1, using both SMC and Calzetti extinction curves. In 
all the samples, the relative accuracy is a^z ^ 0.05 where Az — {zspe — Zphot)/{^ + ^spe)- 

Our combined sample allows to compute the 4400 A rest frame LF in a continuous way 
from z ~ 0.4 up to 2; ~ 3.5. The galaxy luminosity in this rest-frame wavelength is not 
strongly affected by dust, is directly observed up to high redshift in the NIR bands and 
can be easily compared with local estimates of the blue luminosity function. As in paper I, 
we used our /—band selected sample for galaxies at 2; < 1. Indeed, the 4400 A rest frame 
wavelength is within or shorter than the I band, and the LF includes all the galaxies with 
m [4400(1 -I- ^)] < lABilim) of the various samples. This implies that some galaxies from the 
original Iab limited samples are excluded from the LF since they have a red spectrum and 
consequently m[4400(l -|- z)] are fainter than our adopted threshold. We have in this way 
extracted a complete subsample of galaxies having a continuum magnitude threshold at 4400 
A in each redshift interval, independently of their color (assuming galaxies in this redshift 
interval have {R — I)ab > 0). The same procedure has been adopted at higher z using the 
K OT H bands selected sample. In the 1.3 < 2; < 2.5 and 2.5 < z < 3.5 interval, the 4400 
A rest frame wavelength is within or shortward of the H and K bands, respectively. The 
LF includes all the galaxies with m[4400(l -|- z)] = rriAB < ?7T.^s(/im), where ttiab = Hab 
or Kab for 1.3 < 2; < 2.5 and ttiab = Kab for 2.5 < 2; < 3.5. Finally, as in paper I, the 
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rest frame absolute AB magnitude Mb{AB) has been derived from the theoretical best fit 
spectral energy distribution for each galaxy. The same fit provides the photometric redshift 
for each galaxy. In this way the K -correction has been derived for each galaxy from the 
interpolation between the observed magnitudes using the best fit spectrum. Uncertainties 
in this interpolation produce on average small (< 10%) errors in the rest-frame luminosities 
(Ellis 1997; Pozzetti et al. 2003). 

3. The evolution as a function of color and star- formation 

To explore the evolution of the galaxy luminosity function for different spectral types 
we have analyzed first the color distribution of galaxies as a function of redshift. In Fig.l 
we show the distribution of the U — V colors as a function of z for the two samples, the I 
band selected a.t z < 1 and the H, K band selected at2;=1.3 — 3.5. The rest-frame U — V 
color has been used because it includes the 4000 A break and so it is most sensitive to the 
galaxy properties (age, star formation, etc.). Moreover the U — V color is always sampled by 
the multicolor catalog of our composite optical-NIR sample in a large redshift interval up to 
z ~ 3.5. No correction for dust absorption has been apphed to the rest frame U — V color. 

Previous analyses on the luminosity function evolution have adopted a constant rest- 
frame color to separate red/early- type galaxies from the blue/star- forming population (e.g. 
the CFRS sample (Lilly et al. 1995)). However, when extending this criterion to high 
redshifts we can not ignore evolutionary effects in the stellar population of high z galaxies. 
As a first attempt to include this effect we assumed a scenario where the star formation 
history of the red/early galaxy population is controlled by passive stellar evolution from a 
high formation redshift. Thus we adopted the PLE (pure luminosity evolution) evolutionary 
color track of a galaxy whose local color is similar to the average color of a typical SO galaxy 
as a criterion to separate the galaxy populations. In this way we can follow the blueing of 
the early-type population seen at increasing redshifts assuming that only passive evolution 
was active in the past. 

In Fig.l the "SO evolutionary track" is shown as a function of z. A formation redshift 
oi z — 10, an exponential time-scale r = 3 Gyr and solar metallicity have been adopted to 
provide at 2; = the typical color U — V ~ 2.1±0.1 of this class of galaxies (Schweizer 
& Seitzer 1992). The latest version of the standard Bruzual & Chariot (2003) model has 
been adopted. This evolutionary track naturally follows the evolution in redshift of the 
bhic edge of the underdense regions occupied by red objects in the observed bimodal color 
distribution. The color separation increases as the redshift decreases from f/ — K ~ 1 at 
z ~ 2.5 io U — V l.b aX, z 0.7 as shown in Fig.l. Above the cut, the spread in the 
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colors can be represented by passive evolution with short (r = 0.3, 1, 3 Gyr) exponential star 
formation timescales with a high formation redshift {z — 10). 

Although the adopted evolutionary color cut is somewhat arbitrary, we have verified that 
small changes of the parameters do not change significantly the results and conclusions about 
the evolution of the red and blue populations outlined in section 5. The main constraint is 

that the evolutionary color tracks should provide at 2; ^ colors still consistent with that 
of SO galaxies. There is room for example to reduce the formation redshift from 2; = 10 to 
z = 6. The new tracks would produce a blue-shift by about 0.1 in the color cut at each 
redshift (f/ — V" ~ 2 at 2; ~ 0). The resulting luminosity function evolutions of the red and 
blue populations do not change appreciably and the changes in the evolutionary parameters 
of the luminosity functions are well within the errors. 

It is clear however that if we are interested in the analysis of the galaxy evolution as a 
function of the star-formation activity, the use of a color criterion introduces some population 
mixing since it is not possible to distinguish an early-type galaxy from a dusty starburst using 
only the U — V color. 

A further step in this direction can be obtained with a spectral fit to the overall SEDs by 
means of the Bruzual & Chariot (2003) spectral synthesis model. Although some degeneracy 
still remains, we can remove the most obvious starburst galaxies from the locus of early 
type galaxies by considering the specific (i.e. per unit mass) star formation rate m*/m* as a 
function of redshift that is plotted in Fig. 2. The use of m*/m* in this context has a number of 
advantages with respect to using the SFR; for example it is largely independent of cosmology 
and it is less sensitive to the stellar IMF (see e.g. Guzman et al. 1997, Brinchmann and 
Ellis 2000, Brinchmann et al. 2004). The library used to fit the galaxy properties (SFR 
and stellar mass) and to derive photometric redshifts for the faint fraction of the sample is 
described in more detail in paper I and in Fontana et al. (2004). We are aware that the 
absolute values in the m*/m* distribution are subject to uncertainties due for example to 
the estimate of the dust attenuation which depends on the extinction curve adopted and 
to the methods adopted for the mass estimate. We refer in particular to the papers by 
Fontana et al. (2004) and Dickinson et al. (2003) for the detailed description of the methods 
used to derive a reliable estimate of the stellar mass from the SED fitting procedure to our 
multicolor galaxy sample. We have recently quantified a systematic difference between the 
stellar mass derived from the best fit to the SED and that based on the so called "maximal 
age approach" . The masses derived from the best fit method are in general a factor < 2 lower 
(Fontana et al. 2004). Any physical interpretation of the values produced in the observed 
m*/m* distribution should take into account all these uncertainties. 

To explore how much the evolution of the red galaxies can be influenced by the presence 
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of dusty starburst galaxies, especially at high we have separated the population according 
to the evolutionary track of the specific star formation rate m*/m* which corresponds to the 
" SO color track" already adopted. The evolution of the classes separation appears evident in 
Fig. 2. The evolutionary SSFR decreases from log(m*/m*) ~ —9.4 to log(m*/m^,) ~ —10.3 
(in units of yr~^) from z ~ 2.5 to z ~ 0.7. Below the cut, the increasing spread in rn^,/m^, 
with decreasing redshifts can be represented by passive evolution with short (r = 0.3,1,3 
Gyr) exponential star formation timescales with a high formation redshift {z = 10). 

The color/SSFR "SO evolutionary tracks" will be used in section 5 to provide the separa- 
tion of the galaxy populations at different redshifts and to analyze their luminosity function 
evolutions under the assumption that the star formation history of the red galaxies is mainly 
subject to stellar passive evolution. 

4. The bimodal color/stEir- formation distributions as a function of z 

A different selection criterion which is independent of the assumed evolutionary scenario 
is based on the empirical analysis of the color or m*/m* distributions. Previous analyses 
at lower redshifts < 1) by BeU et al. (2004) and Baldry et al. (2004) have shown 
that these distributions are bimodal and so identify two distinct populations, the red/early 
and blue/late galaxies. In the next section we analyse the color and SSFR distributions in 
our composite sample up to 2; ~ 3 to probe the presence of bimodality at high redshifts. 
The bimodality will be used as a further criterion to study the evolution of the luminosity 
functions of the two populations. 

4.1. The bimodal color distributions 

The inclusion of the deep Hubble Deep Fields H and band selected galaxies in our 
composite sample shows that this bimodality continues at least up to 2; ~ 2 — 3 and mainly 
separates blue/star-forming galaxies from the red/early-type population. 

Previous analyses (e.g. Bell et al. 2004) have also shown that there is a weak luminosity 
dependence of the color separation induced by a color-magnitude relation. In Fig.3 this 
relation is shown for the two redshift intervals. To get a first evaluation of this relation we 
have simply fitted the distribution shown in Fig.3 with the sum of two gaussians having 
the mean which is a linear function of the blue absolute magnitude. Each gaussian has a 
constant dispersion and each survey has been weighted with its covering sky area. Since 
the statistics of the red population are poor we have adopted the same dependence on the 
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absolute magnitude for both the blue and red populations. The maximum likelihood clearly 
indentifies the locus of the blue population in the redshift interval z = 0.4 — 1, with a 
mean color which correlates with the absolute magnitude following the relation {U — V) = 
(-0.098 ± 0.03) ■ {Mb + 20) + (1.11 ± 0.02) with a scatter a = 0.24. For the red population 
we find {U - V) ^ (-0.098 ± 0.03) • {Mb + 20) + (1.86 ± 0.03) with a = 0.24. 

In the magnitude range present in each redshift interval, the changes in the average 
colors of the two populations are at most A{U — V) ^ 0.3 from Mb = —17 to Mb = —22. 
In the higher redshift interval 2; = 1.3 — 3.5 the statistics are poor and the significance of the 
correlation is weaker. For the blue galaxies wc derive {U — V) = (—0.062 ± 0.017) ■ {Mb + 
20) + (0.73 ± 0.01) with a = 0.25 and for the red population {U - V) = (-0.062 ± 0.017) • 
{Mb + 20) + (1.65 ± 0.28) with a — 0.42. An intrinsic blueing as a function of redshift is 
apparent for both populations, as expected even in the case of pure passive evolution. At 

Mb = -20 from z ~ 0.7 to z ~ 2.5 it is A{U - V) 0.4, -0.2 for the blue and red 

population, respectively. 

Taking into account the different normalizations of the two gaussian distributions, it is 
possible at this point to find a formal minimum for the sum of the two distributions, which 
allows to separate in an empirical way the blue from the red galaxy populations. This locus 

is represented by the following color-magnitude relation {U — V) = —0.098 ■ (Af^ + 20) + 1.64 
at ^ = 0.4 - 1 and {U - V) = -0.062 ■ {Mb + 20) + 1.31 at ^ = 1.3 - 3.5 (dotted lines in 
Fig. 3). In the low redshift bin our U — V color separation is in good agreement with that 
found in the C0MB017 survey at ^ ~ 0.5 ( (f/ - ~ 1.6 at M^j = -20, after changing h 
and adopting an average Mb — My = 0.5 as derived from the red sample). 

The histograms of the color distribution are shown in Fig.4 in the two redshift intervals. 
The histograms have been projected at the average characteristic magnitude of the two 
Schechter populations in the two redshift intervals for illustrative purpose. The vertical 
lines represent the average separation between the two populations evaluated from the color 
bimodal distribution. 

We note in this context that the high redshift interval where the bimodal distribution 
is measured is rather large as required to keep sufficient statistics. Only few red objects are 
indeed present in the sample ai z > 2 and therefore we can not firmly estabfish the presence 
of a bimodal distribution at these redshifts. However, the color distribution at 2; > 2 (see the 
dotted histogram in the lower panel of Fig.4) is at least consistent with such a bimodality 
being already in place at 2; > 2. 
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4.2. The bimodal rh^/m* distributions 

Considering the corresponding m*/m* distributions, a bimodal shape appears at all red- 
shifts, although its dependence on the rest-frame blue absolute magnitude is only marginally 
significant in the low redshift bin and fully consistent with no correlation in the high redshift 
bin (Fig. 5). Indeed the same fitting procedure with two Gaussians gives (logm,,/m*) = 
(0.036 ± 0.016) • {Mb + 20) - (9.09 ± 0.03) for the "late" type galaxies and (log m*/m*) = 
(0.036 ± 0.016) • {Mb + 20) - (11.5 ± 0.8) for the "early" population in the redshift interval 
z = 0.4 — 1. In the maximum likelihood analysis we have iteratively removed objects above 
3(7 from the resulting Gaussian distributions. We also note that we call "late" and "early" 
the two main populations of the bimodal distribution without any specific reference to their 
morphological properties. 

In the higher redshift interval z = 1.3 — 3.5 we obtain (logm*/m*) = (—0.015 ± 0.040) • 
(Mb + 20)-(8.7±0.1) for the "late" population and (log m*/m*) = (-0.015±0.040) ■ (Mb + 
20) — (11.4±0.2) for the "early" population. A small increase (A(Iog m^/m^) ~ 0.3, 0.1) in the 
mean value of the specific star formation rate with increasing redshift is apparent at Mb = 
—20 for both populations, respectively. This increase of course reflects the corresponding 
blueing of the two populations with increasing redshift, as recently remarked e.g. by Rudnick 
et al. (2003) and Papovich et al. (2004). The errors in log m*/m* derived from the spectral 
fit to each galaxy SED are of the same order (~ 0.15) for almost all the objects in the 
sample and therefore error weighted and unweighted ML analyses produce essentially the 
same result. 

As in the case of the color distribution, it is possible to define a minimum between the 
two gaussian distributions represented by the relation log m*/m* = 0.036 • {Mb + 20) — 10.7 
at z = 0.4 — 1 and log rh^/m* = -0.015- (Mb -|- 20) — 10.6 at z = 1.3 — 3.5, weakly dependent 
on the absolute magnitude. The histograms of the SSFR distribution are shown in Fig. 6 in 
the two redshift intervals as in Fig. 4 for illustrative purpose. 

Thus the much stronger and more significant color-magnitude correlation does not ap- 
pear to be mainly due to appreciable changes of the specific star formation rate with the 
rest frame blue magnitude, but more likely to variations of some other physical parameters. 
Keeping in mind the various degeneracies involved in the derivation of the galaxy physi- 
cal parameters from the color fitting to the synthetic spectral models (see e.g. Papovich, 
Dickinson & Ferguson 2001; Fontana et al. 2004), it is interesting to note that the fitting 
analysis suggests dust reddening as the main responsible for the observed color-magnitude 
correlation in the blue population. Indeed the fraction of the blue galaxies with extinction 
E{B — V)> 0.2 (big symbols in Fig. 3) increases with the blue luminosity. In the redshift 
interval z — 0.4 — 1, considering the blue galaxies with colors within 2a from their average 
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color-magnitude relation, the fraction of dusty galaxies is 36% and 60% for galaxies with 
Mb > —20 and Mb < —20, respectively. In the higher redshift interval 1.3 < ^ < 3.5 the 
fraction increases from 24% at low luminosities to 40% at high luminosity, but the statistics 
is admittedly poor. This is in agreement with previous findings by Shapley et al. (2001) for 
the Lyman break galaxies at 2; ~ 3. 

For the red population the change in the fraction of dusty galaxies from faint to bright 
galaxies is only 5% and it is not sufficient to explain the observed color-magnitude correlation. 
However, the fraction of red galaxies with solar or greater than solar metallicity is higher 
(65%) for high luminosity galaxies {Mb < —20) than for faint galaxies (50%) in agreement 
with what is found in the local early type galaxy sample (see e.g. Kuntschncr et al. 2002). 
This last effect could more likely reproduce the observed correlation even at high redshift. 

Although these results are dependent on several assumptions like e.g. the adopted 
Salpeter IMF, they can help in identifying the main physical parameters that underlie the 
observed relations. Further analysis in this direction is out of the scope of the present work. 

Summarizing the analysis of this section, we have shown with this still limited sample 
that the bimodality in the distribution of the galaxy U — V colors or specific star formation 
rates extends to high redshifts at least up to ;2 ~ 2. A more robust statistical significance 
will be obtained with larger multicolor optical-NIR surveys. 

Despite all the uncertainties, the derived bimodal distribution is clearly representative of 
the presence of two distinct populations and it can be used to make an empirical separation of 
the galaxy populations at different redshifts to analyze their luminosity function evolutions. 
This will be done in the next section. 



5. The evolution of the luminosity functions of field galeixies 

5.1. The statistical analysis 

The luminosity function is computed extending the procedure used in paper I. First of 
all we have applied to our composite sample an extended version of the l/Vmax algorithm 
as defined in Avni and Bahcall (1980) so that several samples can be combined in one 

calculation. Our combination of data from separate fields, with different magnitude limits, 
has been treated computing, for every single object, a set of effective volumes, VmaxU), 
for each j-th field under analysis. For a given redshift interval (^1,^2), these volumes are 
enclosed between Zi and Zup{j), the latter being defined as the minimum between Z2 and 
the maximum redshift at which the object could have been observed within the magnitude 



limit of the j-th field. The galaxy number density 0(M, z) in each {Az, AM) bin can then 
be obtained as: 



«=1 L j ^ 

where is the area in units of steradians corresponding to the field j, n is the number 
of objects in the chosen bin and dV/dz is the comoving volume element per steradian. 

The Poisson error in each LF magnitude bin was computed adopting the recipe by 
Gehrels (1986), valid also for small numbers. The uncertainties in the LF due to the pho- 
tometric z estimates were computed as in paper I, assuming the typical r.m.s. found com- 
paring photometric and spectroscopic z. A set of 50 catalogs of random photometric z was 
produced using the above r.m.s. uncertainties to compute a set of 50 LFs. The derived fluc- 
tuations in each magnitude bin resulted smaller than Poissonian errors and have been added 
in quadrature although the errors are not independent between different bins and could be 
overestimated. Finally, in order to take into account, at least partially, the uncertainties 
due to the field-to-field variance, we have used the maximum between the Poisson error and 
the density variance among the fields in each bin of the luminosity function. The resulting 
field-to-field variance is the main source of error in the LF bin. 

The 1/Vmax estimator for the LF can be affected in principle by small scale galaxy 
clustering although in our case the large redshift bins adopted would reduce the effect of 
any density fluctuation. For this reason a parametric maximum likelihood estimator is also 
adopted which is known to be less biased respect to small scale clustering (see Heyl et al. 
1997). 

The parametric analysis of the galaxy LF has been obtained from the maximum like- 
lihood analysis assuming a Schechter parametric form for the LF. The method used in 
Poll et al. (2003) and in the present paper represents an extension of the standard Sandage, 
Tammann & Yahil 1979 method where several samples can be jointly analyzed and where 
the LF is allowed to vary with redshift. In the case of our composite survey, that is based 
on j magnitude-limited flclds, the maximum likelihood A has been derived from the analysis 
of Marshall et al. (1983) after proper normalization of their probability function: 



N 



A=n 



0(Mi, Zi)dzdM 



1 Ej^ij) ff %dz J^J^-^'^ m, z)dM 



(2) 
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where zi,Z2 are the minimum and maximum redshifts available for each sample, N is 
the total number of objects considered in the composite sample, Mi-^{z) is the absolute 
magnitude value corresponding, at the given redshift z, to the magnitude limit of the j-th 
survey. The value of </>* in Table 1 is then obtained simply imposing on the best fit LF a 
normalization such that the total number of galaxies of the combined sample is reproduced. 
A formal derivation of this equation can also be found in Heyl et al. (1997) where it is also 
shown that the present equation reduces to the standard Sandage et al. formulation in the 
case of a single survey and of a LF which is constant in redshift. 

At variance with paper I we allow the Schechter parameters 0* and M* to vary with 
redshift : 

(j){M,z) = 0.4 ■ ln(lO)0*(z)[lO°-^(^*(^)-^V^"exp[-lO°-^(^*(^)-^)] 

where 0* = 0q ■ (1 + z^ and M*{z) = Mq — 6 ■ log(l + z) equivalent in luminosity to 
L*{z) = Lq ■ (1 + z)'^^'^^. The slope a is kept constant with redshift since it is well constrained 
only by the low- intermediate redshift data (z < 1). 

We used the MINUIT package of the CERN library (James & Roos 1995) to minimize 
—2lnA. Since both the composite sample and the parametric analysis have been improved 
respect to paper I, we show in Fig. 7 and in the first row of Table 1 the results about the 
evolution of the luminosity function of the overall sample. The errors in Table 1 correspond 
to errors for each interesting parameter, independently of the other parameters, obtained 
with a standard threshold of 2AlnA — 1. A detail analysis of the 4-dimensional space is 
out of the scope of the present paper. However, if we consider e.g. the joint 68% confidence 
region for the M*, a or 7, 5 pairs of interesting parameters we obtain, after projection on the 
parameter axes, individual errors that are about 1.6 times larger. 

The inclusion of the deep NIR selected HDF-N sample contributes to a better definition 
of the faint LF slope at intermediate redshifts {z = 1.3 — 2.5). As in paper I, the local 
LFs from the Sloan (Blanton et al. 2001) and 2dF (Norberg et al. 2002) surveys have 
been included for comparison. The difference between the two local LFs show the level of 
uncertainty present in the estimate of the LF even at about the same local redshift interval. 

The evolution of the LF shown in Fig. 7 by the 1/Vmax analysis shows little density evo- 
lution at the faint end with respect to the local values, while at the bright end a brightening 
by about 1 magnitude is apparent. 

If we adopt a Schechter shape of the LF at all redshifts, the maximum likelihood analysis 
shows an average slope, (a ~ —1.3), somewhat flatter than that found in Paper I. Moreover, 
appreciable density and luminosity evolution with redshift arc needed to fit the data, as 
derived from the parameters in Table 1. We note in this respect that the magnitude limits 
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of the samples cut the LF at progressively higher luminosities with increasing redshifts. For 
this reason wc prefer to keep the a slope constant with redshift unless a bad agreement 
with the 1 /Vmax points is found. The mixed luminosity-density shift of the Schechter shape 
provides a luminosity function which keeps an almost constant volume density with redshift 
at the faint end, while at the bright end it provides an appreciable luminosity evolution, in 
agreement with the qualitative results of the XjVmax analysis. 

Fig. 7 also shows the good agreement, better than the la level, between the parametric 
and the l/Knaa; LFs. This implies in particular that the assumed parametric evolutionary 
laws arc an acceptable representation of the binned luminosity function in the overall redshift 
interval. Moreover, the good agreement implies that any small scale clustering does not affect 
appreciably the l/Knax estimate of the LF due to the large redshift intervals adopted. 

This kind of shape and evolution of the luminosity function is the result of the different 
evolutions of distinct galaxy populations. For this reason we perform in the next section 
a separation of the galaxy sample in the blue/late and red/early populations to probe the 
different evolutionary behaviour of the two luminosity functions. 

5.2. Passive evolutionciry tracks £is a criterion to sepetrate red/eeirly and 

blue/late galeixies eis a function of z 

The luminosity functions of the blue and red populations defined on the basis of the 
"SO color track" are shown in Fig.8 in different redshift intervals. The parametric Schechter 
evolution has been derived adopting a mixed luminosity and density evolution with a constant 
slope. The model luminosity function is shown at the average z of the bin for illustrative 
purpose. 

The results of the Maximum Likehhood fits are shown in Table 1. Assuming the 
Schechter shape at all redshifts we found that the average slope a of the two populations is 
clearly different, ranging from a. ~ —1.3, —1.4 for the blue/late galaxies to a ~ —0.5, —0.6 
for the red/early population. This result is already well known from previous works on the 
local and ^ < 1 samples (Lilly et al. 1995, Blanton et al. 2001, Bell et al. 2004). The 
shallower slope for the red/early sample implies that the ratio of the blue over the red pop- 
ulation increases strongly to fainter magnitudes. The slope of the red LF in particular is in 
good agreement with the value of —0.6 ± 0.1 derived in the 2; = — 1 interval by Bell et al. 
(2004) in their analysis of the wider but shallower C0MB017 survey. The 5 and 7 parame- 
ters of the luminosity and density evolution are significantly different from suggesting that 
both evolutions are needed to provide a satisfactory fit. Both populations show comparable 
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luminosity evolution, but the red population is characterized by a stronger density evolution 
with a decline by at least a factor ~ 5 from 2;~0.4tO2;~3. 

A comparison with the shape and evolution of the total LF computed in the previous 
section (see Table 1) shows that both the faint end slope and the amount of density evolution 
of the total LF are mainly driven by that of the blue population. 

To explore how much the evolution of the red galaxies can be influenced by the pres- 
ence in this population of dusty starburst galaxies, especially at high we have separated 
the population according to the evolutionary " SO track" of the specific star formation rate 
m*/m*. 

The resulting luminosity functions are similar to the previous ones with only mild dif- 
ferences (Fig. 9). From the last column of Table 1, where the number of objects in each 
sample is given, it can be noted that about 35% of the red sample with colors redder than 
the " SO color track" shows a spectral energy distribution indicative of a dusty star forming 
population. This fraction is similar to what is found in spectroscopic surveys and does not 
alter the qualitative behaviour about the LF evolutions. 

It is clear that the presence of a mixed luminosity-density evolution in both populations 
implies, in particular for the red/early population, a star formation history more complex 
than provided by the simple PLE model adopted here to separate the two populations, 
where the burst of star formation is located at high redshifts {z > 3). This mixed evolution 
is directly seen in the 1/Vmax representation of the LFs in Fig.8 and Fig.9, which shows 
that luminosity evolution is apparent mainly at the bright end of the LFs, while at the faint 
end the LFs appear affected by density evolution. This behaviour will be discussed more 
quantitatively in section 6. 

5.3. The empirical bimodal color and m*/m* distributions as a criterion to 
sepEirate red/eeirly and blue/late galeixies as a function of z 

In the previous section we have shown that the evolutionary behaviour of the red/early 
population does not follow what expected by PLE models, assuming a passive evolutionary 
scenario for the color or m*/m* separation. A selection criterion which is independent of 
the assumed evolutionary scenario is based on the empirical analysis of the color or m*/m* 
distributions as observed in our sample and already found at lower redshifts by Bell et 
al. (2004) and Baldry et al. (2004). In the previous section we have shown that these 
distributions are bimodal up to 2; < 3 and so identify two distinct populations. The peaks of 
the bimodal color distributions derived in the low and high redshift intervals are dependent on 
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both redshift and luminosity, while the m^jm^ distributions are less dependent on luminosity. 

Using the locus of the minimum in the color and m^jra^ distributions in the 0.4 < z <\ 
and 1.3 < 2; < 3.5 intervals, we have derived the evolution of the luminosity function of 
red/early and blue/late galaxies shown in Fig. 10 and Fig. 11, whose parameters are shown 
in Table 2. In this case, in each of the two redshift intervals, the color separation is a function 
of the rest-frame blue absolute magnitude. At = —20 the color that divides the two 
classes becomes redder by 0.3 mag from z ~ 2.5 to 2; ~ 0.7, similar although a bit smaller 
than predicted by the PLE color track used in the previous section. 

For the red population we have compared our LF with that of Bell et al. (2004) finding 
a general good agreement at relatively bright magnitudes and low redshifts. In particular, 
the slope of the red population is a = —0.76 ± 0.10 and even in this case it is in agreement 
with that found in the COMBO 17 survey. 

We do not consider the effects of uncertainties in the location of the minima of the 
bimodal distributions, but we do not find a change in the qualitative evolutionary behavior 
of the two populations comparing the LF results obtained using the SO evolutionary track 
and the minima in the bimodal distributions. The main quantitative difference consists in a 
stronger density evolution (larger 7) for the red/early population in the bimodal case (since 
we are excluding some objects at high redshifts with respect to the "SO cut"). 

Indeed we can compare the luminosity function of the early population selected to have 
logm*/m* < —10.5,-10.6 nearby the characteristic absolute magnitude M* = —21,-21.5 
(depending on redshift). At these magnitudes we are on the flat part of the luminosity 
function and thus we are less sensitive to changes in volume density due to any luminosity 
evolution. In Fig. 11 we can see that the LF of the "early" sample dechnes by about one 
order of magnitude from 2; = 0.5tO2; = 3. This density evolution is somewhat smaller when 
considering the color-selected red sample because of the presence of some dusty starburst 
galaxies at high redshifts. Indeed ~ 40% of the red sample shows SEDs indicative of a dusty 
star-forming population. 

6. Number and luminosity density evolutions of the red/early and blue/late 

galetxies 

6.1. Observed behaviour 

The values of the parameters of the mixed evolution we have derived for both populations 
are the result of the assumed parameterization of the LF with a Schechter function. To clarify 
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in a less model-dependent way the cosmological evolution of galaxies at different luminosities 
we evaluate directly the volume density of galaxies as a function of redshift in given bins of 
the rest-frame blue magnitude. Wc can also evaluate the associated redshift evolution of the 
luminosity density observed in the rest-frame B band by the same galaxies. 

We show the results both assuming the " SO evolutionary track" and the empirical sep- 
aration derived from the bimodal distribution to separate the two populations. 

Considering the former case we note in Fig. 12 that the cosmological evolution of the 
red/early galaxies appears to depend on their luminosity. Indeed the faint galaxies, typically 
in the range —21.5 < Mb < —20.6 (where Mb ^ —20.6 is at the faintest hmit of the 
population at z ^ 3), decrease by a factor ~ 2 — 3 from z ~ 0.5 to z ^ 2.5. However it is 
interesting to note that the brighter galaxies in the range —23 < Mb < —22 are consistent 
with being constant in the same redshift interval. 

To allow a first simple quantitative comparison with a PLE scenario, wc also show in 
Fig. 12 the expected evolution for galaxies with —23 < Mb < —22 obtained adopting the 
local LF of red galaxies and evolving it in luminosity adopting the "SO evolutionary track" 
shown in Fig.l. The local value derived from the Sloan survey is included as derived from 
Bell et al. (2004). Although the C0MB017 sample has been derived from a slightly different 
color threshold, the intermediate redshift values are in good agreement with our data. 

The PLE behaviour is similar to the observed evolution for 2; < 1, but it overestimates 
by a factor ~ 5 the data at high rcdshifts. This overestimate reduces to a factor ^ 3 if 
compared with the number density derived from the best fit LFs. This would imply that 
dynamical events are important for the brightest red/early population at z > 1. 

The B-band luminosity density derived from the same data is shown in Fig. 13 and has 
the same behaviour of the number density evolution. The data (open circles) represent 
the contribution from red/early galaxies with Mb < —20.6 which are at the faintest limit at 
2; ~ 3. In this way we can explore the evolutionary behaviour over the entire redshift range of 
galaxies within the same luminosity range. At 2; < 1 we have also considered (open triangles) 
a lower limit Mb < —19.5 to compare the luminosity density at intermediate redshift with 
that of Bell et al. (2004) finding a broad consistency within the rather large errors. Our 
values however are about a factor ~ 2 higher since the "SO track" used to separate the 
two populations tends to include in the red population a higher number of galaxies with 
increasing rcdshifts respect to the bimodal class separation used by Bell et al. (2004). The 
same PLE model is also shown in Fig. 13 and the same discrepancy is apparent at 2; > 1. 

Since the evolution is more complex than expected from a PLE scenario, it is interesting 
to use the empirical criterion for the class separation based on the star-formation bimodal 
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distribution. The volume densities and the B-band luminosity densities of galaxies with 
Mb < —20.6 as a function of redshift are shown in Fig. 14. In this case both the early and 
late populations are shown in each panel. 

The cosmological evolutions of the bright "late" and "early" populations appear to have 
a bifurcation at 2; > 1. Indeed while star-forming galaxies slightly increase or keep constant 
their number /luminosity density, "early" galaxies decrease in their number density by about 
one order of magnitude and in luminosity density by a significant factor ~ 5 — 6 (although 
with large statistical uncertainties) from z ~ 0.4 to 2; ~ 2.5 — 3. We have also plotted the 
luminosity density of both populations derived from the best fit LFs extrapolated down to 
the minimum luminosity observed in our low z LF, Mb = —15.8. The evolutionary trend is 
similar to that observed although shifted in normalization by a factor < 2. It is important 
to note that the parametric representation used in the present analysis is valid only in 
the adopted redshift interval 0.4 < 2; < 3.5. Extending the analysis outside this redshift 
interval could require more complex parametric evolutions. In any case, while at 2 < 1 both 
populations give a comparable contribution to the total B-band luminosity density, at z > 1 
the contribution by the red population is much lower. This results in a global luminosity 
density of galaxies with Mb < —20.6 almost constant with redshift in the same z interval. 

Since we have seen in paper I that there is a strong correlation between the observed 
rest-frame blue luminosity and the star-formation rate for the blue/late population, this 
implies that the star-formation activity in star-forming galaxies slightly increases or keeps 
constant up to 2; = 3. On the contrary, the number of galaxies with low specific star - 
formation rate log m*/m* < —10.5 which are representative of an old stellar population with 
high formation redshift z > 3, decreases quickly. 

This behaviour is in qualitative agreement with what expected from a hierarchical clus- 
tering scenario where galaxies with low specific star- format ion are located in the rare peaks 
of high dark matter density at high redshifts. However, a detailed quantitative comparison 
with specific CDM models is required to understand which physical mechanisms govern the 
early evolution of galaxies. A first attempt is done in the next section. 

6.2. Predictions by CDM models 

To understand how the above results fit into current hierarchical models of galaxy forma- 
tion we used our semi-analytic model (SAM). Starting from primordial density perturbations 
in a A-CDM Universe, the model describes the collapse and the merging histories of the DM 
haloes which originate from such perturbations. In previous papers (see, e.g., Menci et al. 
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2002) we have described our model in detail; it includes the dynamical friction and the bi- 
nary aggregations, determining the fate of clumps included into larger DM haloes along the 
hierarchy of mergers. The baryonic processes (gas cooling, star formation, Supernovae feed- 
back) are linked to the dynamical history of the galaxies through the simple recipes adopted 
by current SAMs (Kauffmann et al. 1993; Cole et al. 1994; Somerville & Primack 1999; Poll 
et al. 1999; Wu, Fabian & Nulscn 2000; Cole et al. 2000). 

In this paper we are interested in the whole distribution of colors for any given galactic 
mass and redshift. In turn, such distribution depends on the different galaxy merging histo- 
ries leading to a given galactic mass. Thus, wc updated our model to describe the scatter in 
the galactic properties (star formation rate, colors) due to different merging histories, which 
could not be treated with our earlier model (Menci et al. 2004). To this aim, we built a 
Monte Carlo version of our SAM, where different histories are generated by extracting the 
merging probability according to the prescription of the extended Press- Schechter (EPS) the- 
ory; the corresponding galactic mass function is computed from such numerical realization 
(at variance with the method used in Menci et al. 2002; 2004). 

The properties of the gas and stars contained in the galactic DM clumps are computed 
as in Menci et al. (2004). The only differences are that now i) stars are allowed to form 
with rate = rric/t^, with the star formation timescale assumed to be proportional to 
the disk dynamical time = qrd/vd through the proportionality constant q, which is left 
as a free parameter; ii) the mass Arrifi returned from the cool to the hot phase is estimated 
as Arrih = Esn/v'^sc where Esn = 10^^ ergs eg is computed from the mass Am^, 

of stars formed assuming ?7o ~ (3 — 5) ■ 10^^ as the number of SNae per unit solar mass 
(depending on the assumed IMF), and eo = 0.01 — 0.5 as the efficiency of the energy transfer 
to the interstellar medium whose value is highly uncertain (see, e.g, Kauffmann, White & 
Guiderdoni 1993; Somerville & Primack 1999) and it is considered as a free parameter. The 
model free parameters q and eo are chosen as to match the local B-band galaxy and the 
TuUy-Fisher relation. We adopt a Salpeter IMF with a standard value eo = 0.1 

At each merging event, the masses of the different baryonic phases are replenished by 
those in the merging partner; the further increments Arric, Am*, Am^ from cooling, star 
formation and feedback are recomputed on iterating the procedure described above. In our 
MonteCarlo code we adopt the description of starbursts triggered by major/minor merging 
and galactic fly-by following the formulation in Menci et al. (2004). 

At higher z, the model results match the observed K-band LF (Pozzetti et al. 2002) up 

to 2; = 1.5 and the UV LFs measured for Lyman-break galaxies at 2; = 3 — 4. The global 
star formation history predicted by the model is consistent with that observed up to redshift 
^ a; 4 — 5 (Fontana et al. 2003a; Giavalisco et al. 2004), and the predicted stellar mass 
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functions fit the observed ones (Cole et al. 2001; Fontana et al. 2004) up to 2; = 1.5. 

Given the above agreement concerning the average properties of the galaxy population 
(i.e., the properties averaged over all the merging histories) the comparison of the SAM 
results with the LF of populations with different m*/m* constitutes a higher-order probe of 
the formation history of galaxies which allows us to test also the properties of populations 
with particular merging histories. 

The predictions by our CDM model concerning number and luminosity densities are 
shown in Fig. 14 for both populations (late and early) adopting the same empirical separation 
derived from the observed bimodal distribution. Indeed the model shows a clear bimodal 
distribution pointing toward the same color dichotomy found in the data. A detailed analysis 
of the model bimodal distribution will be presented in a subsequent paper. It is to note that 
since the red/early luminosity functions at all redshifts are rather fiat around Mb ~ —20.5, 
the comparison with the model predictions depends only weakly on the magnitude cut. 

We have to keep in mind that the threshold in the m^,/m=^ adopted to separate the 
two populations could be subject to systematic underestimates especially concerning the 
stellar mass derived from the observed SEDs, as discussed in sect. 3; in particular, the 
"maximal age approach" yields stellar masses a factor 2 larger (Fontana et al. 2004); thus 
we conservatively adopt in the predicted distribution of log{m^,/m^) a spread ranging from 
the default threshold (-10.5, corresponding to the best fit observed stellar mass estimate) 
up to an offset of 0.3. The assumed uncertainty in the threshold propagates to the model 
predictions concerning the number/luminosity densities, as represented by the shaded regions 
in fig. 14. In principle, an additional source of uncertainty is represented by modelling the 
dust extinction curve, which affects the derived SFR. However, for galaxies with m*/m* 
around the threshold value, the low amount of dust makes this uncertainty much lower than 
that related to the stellar mass. 

Figure 14 shows a broad agreement between the observed and the predicted evolutions 
of the early and late type galaxies with Mb < —20.6. In particular, a decrease in the 
number/luminosity density of red/early galaxies by a factor 3 — 4 respect to the local value 
is predicted at 2 < 2; < 3. Thus, finding red objects at high z is not inconsistent — by itself — 
with hierarchical clustering scenarios. This is because the clumps which end up in massive, 
bright objects are predicted to have older stellar populations (and hence redder colors) on 
average since they formed in rare, high density (biased) regions at high redshifts. The higher 
densities present in these regions have provided enhanced star formation, rapidly exhausting 
the cold gas reservoir. 

Thus, at the level of details of the present SAMs it is remarkable that a clear distinction 
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between the evolutionary trends of the two populations is predicted by this class of models. A 
detailed study of the link between the different merging histories leading to the two different 
populations will be presented in a future paper. 

A similar behaviour is also expected in the hydrodynamic simulations by Nagamine 
et al. (2004) where an appreciable number of massive (and hence mostly red) galaxies is 
predicted at z ~ 3. 

In this comparison we should however keep in mind that several uncertainties affect 
also the data points. Field to field variance due to clustering could affect in particular 
the number of early type galaxies observed at high redshifts in small fields. Moreover, the 
average fraction of early type galaxies in the red population (~ 60%; sec last column of 
Table 2) is statistical in nature and based on the fitting of the SED distribution with the 
spectral synthesis models. However, unless this fraction results much different after detailed 
spectroscopic analysis of the red sample, the conclusion is that the present CDM models 
are able to predict the presence of a number of early type galaxies up to 2; = 3.5 which is 
comparable to that observed. 



7. Summciry 

We have used a composite sample of /, H, K selected galaxies to derive the evolution 
of the B-band luminosity function of galaxies separated by their rest-frame U — V color or 
specific (i.e. per unit mass) star-formation rate up to 2; ~ 3. 

The observed color or SSFR distributions appear bimodal up to 2; ~ 2.5 — 3, extending 
to high redshift the same behaviour recently found at low and intermediate redshifts. In 
particular, the bimodality in the SSFR distribution depends weakly on the luminosity of 
the sources in the luminosity interval probed by the present sample and evolves little from 
z ~ 0.5 to 2 ~ 3. 

We adopted two different criteria to classify the galaxies and to study the evolution of 
the luminosity function of the red/early and blue/latc populations. 

Focusing on the evolution of the early galaxies, we first adopted the null hypothesis that 
bright red/early galaxies are mainly affected by substantial passive luminosity evolution. For 
this reason we have adopted the color/SSFR passive evolutionary tracks of a galaxy whose 
local color is similar to that of a typical SO galaxy, to separate the two populations as a 
function of redshift. The resulting qualitative behaviour of the luminosity function does not 
depend significantly on the parameter used for the selection, U — V color or SSFR, implying 
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that the contamination of dusty star-forming galaxies in the red sample, although substantial 
in number, does not modify appreciably the global evolutionary properties of the population. 

Adopting a Schechter shape, the LF slope of the blue/late population is ~ —1.4 up 
to 2; ~ 2.5, while that of the red/early population is flatter (~ —0.5) extending to higher 
redshifts the similar difference found at 2; < 1. The evolution of the luminosity function of 
blue and red galaxies is described in parametric form by the presence of both luminosity 
and density evolution. The evolution of the blue population is mainly driven by luminosity 
evolution, while appreciable density evolution characterizes the red/early population. In 
particular, a decline by a factor of the order of 5 results from z ~ 0.4 to 2; ~ 3. This 
implies that the red/early population as a whole can not be simply described by simple pure 
luminosity evolution models, but major dynamical events could act at high redshifts. 

For this reason we also adopted an empirical criterion to separate the two galaxy popula- 
tions, based on the presence of a minimum in the bimodal observed color/SSFR distribution 
in each redshift interval. The evolutionary difference between the two populations is even 
more marked in this case, with a density decline of the early population by about one order 
of magnitude from 2; ~ 0.4 to 2; ~ 3 at Mb ~ —21, —21.5. 

The amount of mixed luminosity-density evolution depends on the assumed parame- 
terization of the LF. We have also evaluated directly the volume and luminosity density of 
galaxies as a function of redshift in given intervals of absolute magnitudes. We have adopted 
both selection criteria to separate the two galaxy populations. Assuming the "SO passive 
evolutionary cut" we have found that while the faint red/early population slightly decreases 
with redshift, the brightest fraction with Mb < —22 is consistent with being constant with 
redshift up to 2; ~ 3. A simple PLE model where galaxies are formed at high redshifts 
overestimates by a factor of the order of 5 the bright fraction at 2; ~ 3. 

Adopting the empirical criterion we have found that the volume and luminosity densities 
of the early and late populations having Mb < —20.6, which is the faint limit in the highest 
redshift bin, appear to have a bifurcation at 2; > 1. In particular, while star-forming galaxies 
slightly increase or keep constant their luminosity density, "early" galaxies decrease in their 
luminosity density by a factor ~ 5 — 6 from 2;~0.4to2;~2.5 — 3. 

A comparison with our latest version of hierarchical CDM models shows that a clear 
distinction between the different evolutionary histories of the early and late type galaxies is 
predicted, in good agreement with the observations where the contribution in the luminosity 
density of the early population is 10 — 20 times below that of the late population at 2; ~ 3. It 
is clear however that a more detailed comparison requires larger surveys to reduce the effects 
of the field-to-field variance and at the same time a better description in the CDM models 
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of the physical mechanisms governing the different cosmic histories of the two populations. 

We thank the anonymous referee for her/his helpful comments. We also thank A. 
Cimatti for the use of the K20 survey data. 
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Table 1: Luminosity function parameters 





0S 


M* 


a 


S 


7 


N 


Total 


0.0056 


-20.93 ±0.23 


-1.29 ±0.04 


1.48 ±0.52 


-0.65 ±0.19 


1434 


<U- V{SO) 


0.0042 


-20.31 ±0.26 


-1.38 ±0.05 


2.35 ±0.58 


-0.52 ±0.23 


1142 


>U- V{SO) 


0.0106 


-20.23 ±0.28 


-0.46 ±0.10 


2.72 ±0.66 


-2.23 ±0.28 


292 


> fn/m^SO) 


0.0050 


-20.55 ±0.25 


-1.33 ±0.05 


1.96 ±0.56 


-0.60 ±0.22 


1242 


< 7h/m{S0) 


0.0044 


-20.66 ±0.36 


-0.63 ±0.11 


2.08 ±0.82 


-1.80 ±0.34 


192 



Table 2: Luminosity function parameters 
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< U — V{bimodal) 


0.0040 


-20.47 ±0.26 


-1.36 ±0.05 


2.30 ±0.57 


-0.45 ±0.22 


1173 


> U — V{himodal) 


0.0110 


-20.49 ±0.38 


-0.76 ±0.10 


2.70 ± 1.00 


-3.10 ±0.38 


261 


> fn / m{himodal) 


0.0050 


-20.51 ±0.24 


-1.32 ±0.05 


2.15 ±0.53 


-0.56 ±0.20 


1282 


< m/m{himodal) 


0.0110 


-20.23 ±0.41 


-0.53 ±0.13 


3.49 ± 1.01 


-3.65 ±0.47 


152 
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Fig. 1. — Rcst-framc U — V a^s a function of redshift. The continuous curve represents the 
color evohition due to the passive evolution of a galaxy formed at 2; = 10 and with a star- 
formation time-scale r = 3 Gyr. The Bruzual & Chariot (2003) code has been adopted. The 
resulting color at 2; = is typical of an SO galaxy. Dotted and dashed curves are derived 
adopting t — 1 Gyr and r = 0.3 Gyr, respectively, and could represent the reddest fraction 
of the sample. 
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Fig. 2. — Specific star formation rate irn/m as a function of redshift. The curves are obtained 
from the same models as in Fig. 1 
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Fig. 3. — Color magnitude relation in two redshift intervals. Black filled dots are galaxies 
with E{B — V") > 0.2. Small dots are galaxies with E{B — V) < 0.2. Thick lines represent 
the best fit relations for the red and blue populations. Dotted fines represent the locus of 
the minimum between the two distributions used for the class separation. 
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Fig. 4. — Color histograms in two redshift intervals. The continuous vertical lines represent 
the average separation evaluated from the color bimodal distribution. The lines are evaluated 
at the average redshift and at the average characteristic magnitude of the two Schechter 
populations (M* ~ -19.9 at 0.4 < z < 1 and M* ~ -21.2 at l.S < z < 3.5) for illustrative 
purpose. The dotted histogram represents the color distribution of galaxies in the redshift 
interval 2 < z < 3.5. 
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Fig. 5. — Specific star formation rate m/m as a function of Mb at different redshift intervals. 
Thick and dotted lines as in Fig. 3 
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Fig. 6. — SSFR histograms. The continuous vertical lines represent the average separation 
evaluated from the SSFR bimodal distribution. The lines are evaluated at the average 
redshift and at the average characteristic magnitude of the two Schechter populations (M* ~ 
-19.7 at 0.4 < 2; < 1 and M* ~ -21.2 at 1.3 < 2; < 3.5) for illustrative purpose. 
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Fig. 7. — Total luminosity function as a function of rcdshift. Continuous curves come from 
the Maximum Likelihood analysis. The local LFs from the Sloan and 2dF surveys are also 
shown as short and long-dashed curves, respectively, as in paper I. 
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Fig. 8. — Luminosity function of blue (left) and red (right) galaxies separated in the rest- 
frame U — V color following the " SO evolutionary track " (see text for details) . Continuous 
curves come from the Maximum Likelihood analysis. The luminosity function in the lowest 
redshift bin 0.4 < 2; < 0.7 is also shown for comparison in all panels (dotted curve). 
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Fig. 9. — Luminosity function of "late" (left) and "early" (right) galaxies separated on the 
basis of their specific m*/m* star formation rate following the "SO evolutionary track". The 
luminosity function in the lowest redshift bin 0.4<2;<0.7is also shown for comparison in 
all panels (dotted curve) . 
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Fig. 10. — Luminosity function of "late" (left) and "early" (right) galaxies separated follow- 
ing their rest-frame U — V colors as derived from the minimum in their bimodal gaussian 
distributions. The luminosity function in the lowest redshift bin 0.4 < z < 0.7 is also shown 
for comparison in all panels (dotted curve). 
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Fig. 11. — Luminosity function of "late" (left) and "early" (right) galaxies separated follow- 
ing their specific m*/m* star formation rate as derived from the minimum in their bimodal 
gaussian distributions. The luminosity function in the lowest redshift bin 0.4 < z < 0.7 is 
also shown for comparison in all panels (dotted curve). 
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Fig. 12. — Density evolution as a function of redshift for red/early type objects separated on 
the basis of the "SO evolutionary tracks". Circles are galaxies in the range —23 < Mb < —22 
and in the same redshift intervals used to compute the LFs. Open triangles are galaxies in 
the range —21.5 < Mb < —20.6. Filled triangles come from the Schechter LFs of Bell et 
al. (2004). The thin continuous curve comes from the PLE model discussed in the text and 
computed at Mb — —22.5. The thick continuous curve comes from the integration of the 
best fit LFs in the range —23 < Mb < —22. Upper panel for the color selected galaxies, 
lower panel for galaxies selected according to their specific star formation rate, as in Figs. 7,8. 
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Fig. 13. — Luminosity density in the rest-frame B-band of the red/early type galaxies selected 
as in Fig. 12. Galaxies separated by their rest-frame U — V color are shown in the upper 
panel, while those separated by their rh^/m^ are shown in the lower panel. The open circles 
represent the contribution from bright galaxies with Mb < —20.6 which are at the faintest 
limit of our sample at z ~ 3. Open triangles represent the contribution from galaxies with 
Mb < —19.5 for comparison with the Bell et al. 2004 data (filled triangles) at 2; < 1. The 
curve represents the passive evolution model shown in the previous figure for galaxies with 
Mb < -20.6. 
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Fig. 14. — Number and B-band luminosity densities of early (filled squares) and late (open 
squares) type galaxies separated by their bimodal m^/m^ distribution. Only galaxies with 
—24 < Mb < —20.6 are considered. The B-band luminosity density of the bright early 
plus late type galaxies is shown as small filled circles in the lower panel. The continuous 
lines represent the "total" luminosity densities of the early and late type galaxies obtained 
after integration of the best fit LFs down to Mb — —15.8. The shaded regions represent 
the uncertainties of the predictions of our CDM model due to the spread in the m*/m* 
threshold adopted and to the dust extinction curves (Small Magellanic Clouds, Milky Way 
and Calzetti) affecting the luminosity cut. 



